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Abstract 

The number partitioning problem can be interpreted physically in terms of a 
thermally isolated non-interacting Bose gas trapped in a one-dimensional harmonic 
oscillator potential. We exploit this analogy to characterize, by means of a detour to 
the Bose gas within the canonical ensemble, the probability distribution for finding 
a specified number of summands in a randomly chosen partition of an integer n. It 
is shown that this distribution approaches its asymptotics only for n > 10 10 . 

Consider the decompositions of a natural number n into natural summands, without 
regard to order. Let $(n,M) denote the number of such partitions which consist of 
M parts, and Q(n) = J2m=i $( n > M) the total number of partitions. For n = 4, for 
instance, we have 

4 = 1+1+1+1 
= 2+1+1 

= 2 + 2 = 3 + 1, (1) 

hence $(4,4) = 1, $(4,3) = 1, $(4,2) = 2, $(4,1) = 1, adding up to 0(4) = 5. It 
is known that f2(n) grows exponentially with yfn Q, so that the enumeration of the 
individual partitions soon becomes impractical when n gets larger. It is then useful to 
focus on the distribution 

$(n,M) , 

Pmc (n, M) = K > ' (0 < M < n) , (2) 

which gives the probability for finding M summands in a randomly chosen partition 
of n. For moderately large n, this distribution can be computed numerically with the 
help of the recursion relation 

min{n-M,M} 

$(n,M)= $(n-M,fc); (3) 

fc=i 

fig. H depicts the results for n = 1000 and n = 5000. 
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Figure 1: Exact "micro-canonical" probability distribution (0) for finding M summands 
in a randomly chosen partition of n, for n = 1000 (dashed) and n = 5000 (full line). 

The number partitioning problem [Q, ^, |4|, |j finds profound applications in various 
areas of statistical physics, ranging from lattice animals || 0] over combinatorial opti- 
mization H to Fermion-Boson transmutation Q. Therefore, it is of substantial interest 
to characterize the distribution (||) for asymptotically large n: Does it, for instance, 
become Gaussian? 

In this Letter we tackle this number-theoretical question within a physical frame- 
work. There is a one-to-one correspondence between the individual partitions of an 
integer and the individual microstates of a gas of ideal Bose particles stored in a one- 
dimensional harmonic oscillator potential with frequency ujq. If the total excitation 
energy E of the gas amounts to n oscillator quanta, E = uHloq, then each partition 
of n labels one possibility for distributing E among the particles. For n = 4, the first 
line in eq. (Q) indicates a microstate with four quanta %ujq bestowed on four different 
particles, the second line indicates another microstate where one particle carries two 
quanta while two other particles account for the remaining ones, and so forth. For any 
value of E, we assume that the number of particles be at least as large as the number of 
quanta, so that no restriction (with respect to the number of summands) on the parti- 
tions occurs. Hence, each excited Bose particle gives a nonzero summand in a partition 
of E/(Hujo) = n, while the remaining ground-state particles correspond to additional 
zeroes. 

Thus, the partitioning problem is mapped to microcanonical statistics: Given the 
total energy (the number to be partitioned) of the thermally isolated gas, the task is to 
count all accessible microstates. The physical model now suggests to consider first the 
simpler canonical version of this problem, i.e., a gas of infinitely many, harmonically 



trapped ideal Bosons in thermal contact with a reservoir of temperature T [ 10 1 . Then 
the microcanonical distribution (||) is replaced by its canonical counterpart 

Pen [pi M) Ene -6n 0(n) > W 

where b = fiujQ/{k-QT) quantifies the inverse temperature, made dimensionless with the 
quantum hujQ and Boltzmann's constant ks- Within the canonical ensemble, the analysis 
starts from the M-particle partition functions 

oo 

Z M (b) = w(n, M) exp(-M6/2 - bn) , (5) 

n=0 
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where the weight u(n,M) is the number of possibilities for distributing n quanta over 
up to M Bosons. Since $>(n,M) counts the number of possibilities for distributing the 
n quanta over exactly M particles, we have 

u(n,M) -u(n,M- 1) = $(n,M) . (6) 



Following textbook practice we proceed from the canonical to the grand canonical 
ensemble by introducing the fugacity z, and defining the grand partition function 

oo 

E(b,z) = £ (W>/ 2 ) M Z M (6) 

M=0 

oo oo 

= £ * M w ( n > M ) exp(-6n) 

M=0 n=0 
oo ^ 

n 1-zexpf-M ' (7) 



u=0 



If we now multiply this function ([?]) by (1 — z), so that the ground-state factor (u = 0) 
is removed, eq. @ provides the necessary link to the desired quantities 3>(n, M): 

E cx (b,z) = (l-z)E(P,z) 

oo oo 

= zM $ ( n ' M ) exp(-ftn) 

M=0 n=0 
oo ^ 

Rl-zexpi-bv) • (8) 

Thus, the grand partition function H ex (6, z) of an ideal Bose gas with amputated ground 
state generates the microcanonical weights <3?(n,M). According to probability the- 
ory [12], the logarithm of H ex (6, z) then generates the cumulants K^{b) of the canonical 
distribution (||): 

lnS ex (M = £^(ln,r • W 

z/=0 ^ 

The first cumulant Ken (£>) is the expectation value of the number of excited particles 

(2) 

at the given temperature, the second cumulant k ci1 (6) its mean-square fluctuation; in 
general, K C n is related to the k-th. central moment of the underlying probability distri- 
bution [12|. In particular, nin = for k > 3 if that distribution is Gaussian. 

It is crucial that these canonical cumulants can easily be calculated in the rele- 
vant temperature regime k^T S> Uojo, that is, for b <C 1: Starting from the product 
representation (g), one derives [ 13 1 the exact formula 

1 rr+ioo 

«$(&) = 7T- / d * r W C(* + 1 ~ *) , (10) 

Z7T7 Jr—ioc 

where T(t) denotes the Gamma function, and is Riemann's Zeta function. This 
formula ( |l0|) allows one to determine the asymptotic expansion of the cumulants from 
the residues of the integrand. Doing the math yields 

~ K ln l + 7 ) + I-li4 +0(63) (12) 
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7T~ 

2C(3) 

7T 4 

156 4 



1 1 

2b + 24 

1 

~ 240 ' 



(13) 
(14) 
(15) 



where 7 ~ 0.57722 is Euler's constant; all higher cumulants are obtained in the same 
manner. 

For applications to the partitioning problem, however, we have to abandon the notion 
of an externally imposed temperature, and to return to a thermally isolated gas. To 
this end, we write the generating function (||) as 



E 

v=0 



e~ bv Y(u, z) 



(16) 



where the series 

Y{v,z) = £ z M ®{v,M) (17) 

Af=0 

is of central importance, since the microcanonical weights $(1/, M) directly figure as 
coefficients. Hence, its logarithm generates the cumulants Kmc(^) of the microcanonical 
distribution (|2|). This function Y(n,z) describes an ideal Bose gas which exchanges 
particles, but no energy with a reservoir, and thus coincides with the partition function 
for the recently introduced Maxwell's Demon ensemble jE3|. Writing e~ b = x, we extract 
Y(n,z) from the series (16) by means of a complex contour integral, 

E ex (b(x),z) 



Y(n,z) 



1 

27d 



dx ■ 



(18) 



where the path of integration encircles the origin of the complex x-plane counter- 
clockwise, and evaluate this integral within the usual saddle-point approximation 
The saddle point bo (z) is determined by setting the logarithmic derivative of the inte- 
grand to zero, resulting in the equation which links energy with temperature, 



n + 1 



^lnH ex (M 



60 (*) 



Within the Gaussian approximation, one is then led to 



lnY(n,z) = In E cx (b (z), z) + nb (z) 



In 2vr 



1 



In 



d_ 

db 



In H ex (b, z) 



(19) 



(20) 



60 0) 



from which the desired microcanonical cumulants are obtained by further differentiation, 



(n) = [ z— ) In Y(n, z) 



(21) 



2=1 



The calculations now are straightforward, but tedious, because of the z-dependence of 
the saddle point. We omit the technical details [16| and state the result for k = 1: 
When the thermally isolated gas carries the excitation energy E = nhojQ, with n 3> 1, 
the expectation value of the number of excited Bose particles takes the form 



,(1) 



lD 2 K&Hh) 



+ 



D 2 Kin(bi) 



1 + 



2 D 2 «£ ) (6i) 



(22) 
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where kS' (b) = lnH ex (fc, 1), D denotes the derivative with respect to 6, and b\ = 6o(l) 
has to be taken as function of n. This latter task is achieved by inverting the saddle-point 
equation (|i~9|) for z = 1, yielding 



1 



'6n 



7T 



2tt 2 



+ 0(n- 1 / 2 ) 



(23) 



Hence, the canonical expectation value (12), expressed in terms of the scaled energy n, 
reads 



7T 



In 



6ra 



7T 



+ 7 



+ 



2vr 2 



In 



6n 



7T 



+ 7 + 1 + 



7T 



6 



+ 0{n- 1 / 2 ) . (24) 



The difference between the canonical and the microcanonical expectation value then 
follows from eq. (p2[), utilizing the explicit expressions ( |TTj ) and ( |i"2| ) of the canonical 
cumulants k[? and Ken : 



In 



6n 



7T 



+ 7 



0(n 



(25) 



Thus, we finally obtain the desired asymptotic formula for the expectation value of the 
number of summands in a randomly chosen partition of a large integer n: 



7T 



In 



6n 



7T 



+ 7 



+ 



2vr 2 



2 In 



6ra 



7T- 



7T 



+ 2 7 + 1 + — 

D 



+ 0(n^ x / 2 ) . (26) 



For example, for n = 1000 this expression © yields k^(IOOO) = 94.8073..., while the 
exact value is 94.82177..., so that the error is only about 0.015%. 

The above calculation of the expectation value kS(ii) illustrates the general strategy 
for computing an arbitrary microcanonical cumulant from eq. (|2~I|): Starting from the 
saddle-point approximation (]20| ) to the generating function In Y(n, z), the k-th cumulant 
Kml(n) is expressed in terms of derivatives of canonical cumulants Ken (b), with < £ < k, 
which, in their turn, are obtained explicitly from the integral fromula (|i~0|). In the same 
manner, the r.m.s. -fluctuation of the number of summands is determined as 



o(n) = ( K £l(n)) 



1/2 



n 



3VQ 
2^3" 



In 



6n 



7T 



+ 7 + 1 



+ 0{n- 1 ' 2 ). 



(27) 



Of particular interest are the coefficient 71 (n) of skewness, and the coefficient 72(H) of 
excess (or kurtosis) [O], 



71 (n) 



Kmc (l) 



3/2 



and 72 (n) 



2 ' 



which quantify the deviation of the number partitioning distribution @ from a Gaussian; 
in the Gaussian case, both 71 and 72 are equal to zero. Determining the third and fourth 
cumulant as outlined above, we find 



71 (n) 



+ 



1.1395 + 
1 r 



1 



0.10128 [ln(n)] z - 0.37376 ln(ra) - 1.7078 

0.0075008 [ln(n)] 4 + 0.025681 [ln(n)] 3 
+0.020024 [ln(n)] 2 - 0.23028 ln(n) - 0.56984] + C(n" 3/2 ) (28) 



71 



5 



1.20 



1.15 



1.10 




(a) 



2.55 



2.45 



4 6 

log 10 (n) 



10 



2.35 




4 6 8 

iog 10 (n) 



Figure 2: (a) Skewness 71 (re) and (b) excess 72(re) of the number partitioning distribu- 
tion (||). The solid lines indicate numerically computed, exact data; the short-dashed 
lines are the predictions of the asymptotic formulae (f28|) and (|29|), respectively. The 
horizontal lines mark the limiting values. 



and, again including terms of order 0(n 1 

1 



72 (re) 



2.4 + 



1 

+ - 

n 



0.28440 [ln(n)r - 0.56714 ln(re) - 10.064 

0.025276 [ln(n)] 4 + 0.022329 [ln(re)] 3 
0.33809 [ln(n)] 2 + 0.73538 ln(re) + 3.78631 + 0(n" 3/2 ) 



(29) 



Thus, 7i(n) and 72 (re) approach nonzero constants, so that the distribution (|2|) remains 
non-Gaussian: 



lim 71 (re) 



12^6C(3) 



7T' 



1.1395 , 



lim 72 (n) = — . 

n— >oo 5 



However, due to the vexating logarithmic corrections in eqs. (|28| ) and (|29j), this asymp- 
totic behaviour is still masked for merely moderately large n: Figure ^ (a) depicts exact 
values of 71 (n), computed numerically with the help of eq. (§), together with the pre- 
diction of the asymptotic formula (28); fig. |2] (b) shows the same comparison for 72 (re). 
It should be noted that the exact evaluation of the recursion relation @ requires a 
substantial amount of computer memory and therefore becomes quite demanding when 
re is of the order of 10 5 , say, while the limiting values of skewness and excess are well 
approached only for n > 10 10 . 

To conclude: While the number partitioning problem is essentially microcanonical 
in nature, so that one associates "temperature" to natural numbers n on the basis of 
their entropy lnO(n), the equivalent problem of harmonically trapped ideal Bosons is 
approached exactly within the canonical ensemble, when temperature is imposed by an 
external heat bath. By means of such a detour to the canonical ensemble, we have 
characterized the number partitioning distribution in terms of its coefficients of skew- 
ness and excess. Central to our approach is the fact that statistical mechanics concepts, 
such as the partition function, have an intrinsic number-theoretical meaning j2|. Even 
subtle differences between the two ensembles come into play here, such as the usu- 
ally neglected difference between microcanonical and canonical expectation values; see 
eqs. (22) and (|25|). Our results show that the number partitioning distribution adopts 
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its asymptotic shape only for n > 10 10 , so that numerical simulations which inherently 
rely on partitions might not reach the proper asymptotics. The analytical method we 
have employed for computing microcanonical cumulants can be generalized to Bosons 
stored in different types of traps, and thus allows one to study the statistical mechanics 
of thermally isolated Bose gases. 

We thank M. Block, K. T. Kapale, G. Schmieder, and M. O. Scully for useful dis- 
cussions. 
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